A simple encoding of a quantum circuit amplitude as a matrix permanent 
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A simple construction is presented which allows computing the transition amplitude of a quantum 
circuit to be encoded as computing the permanent of a matrix which is of size proportional to the 
number of quantum gates in the circuit. This opens up some interesting classical monte-carlo 
algorithms for approximating quantum circuits. 

PACS numbers: 03.67.-a, 03.67.Lx, 73.43.Nq 



In a recent article pQ Loebl and Moffatt gave a method 
for expressing the computation of the Jones polynomial 
of a braid in terms of a matrix permanent. Although 
computing permanents is believed difficult (#P-complete 
in the language of complexity theory), there exist proba- 
bilistic algorithms [2 which sample the permanent, and 
this suggests some interesting new classical algorithms 
for estimating the output amplitudes of quantum circuits 
because evaluating the Jones polynomial at certain roots 
of unity is BQP-complete [3 j. The route to encoding a 
quantum circuit as the Jones polynomial of a knot, and 
then as a matrix permanent, is somewhat complicated - 
the purpose of this article is to present a simpler con- 
struction. 

We restrict to quantum circuits built from Toffoli and 
Hadamard gates, which are universal [4]. We rely heav- 
ily on the construction of Dawson et al [5]. There is 
it shown how the transition amplitude for such a quan- 
tum circuit is equivalent to counting the number of so- 
lutions of a GF(2) (i.e. XOR-AND) polynomial over 
some binary valued variables. More precisely, the re- 
sults of [5] imply the following: Given a quantum circuit 
U and in/output computational basis states |in), |out) 
the amplitude (out|i7|in) can be expressed as the differ- 
ence in the number of solutions to a GF(2) polynomial 
over (roughly) as many boolean variables as there are 
Hadamard gates in the circuit. It is perhaps easiest to 
explain the construction using an example, such as in 
Figure 1. The a^,^,.. etc are boolean variables, which 
we imagine travelling along the qubit lines. Every time 
the qubit goes through a Hadamard gate we create a new 
such variable, and whenever a variable Z{ travels through 
the target of a Toffoli gate we replace it by z^xiyi where 
Xi,yi are the variables at the control lines of the Toffoli 
gate, as indicated. 

Having labelled the circuit with these variables, we 
then create the function f(x) by taking the sum (mod 
2) of the product of every pair of variables on either side 
of a Hadamard gate. For the example of Fig. [I] we obtain: 
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FIG. 1: Mapping from a standard Toffoli- Hadamard circuit 
to counting solutions of a GF(2) polynomial 



If we are interested in, for example, the amplitude 
(0011 1 U 1 0000) we then fix the input and output variables 
of / accordingly: in this case we would set a\ = bi = c\ = 
d\ = <i4 = 65 = 0, C5 = d 3 = 1 and / simplifies to: 

f(x) = a 2 a 3 © b 2 b 3 © b 3 b A © d 2 c 4 6 4 © c 2 c 3 © b 2 a 2 c 3 © 
C3C4 © c 4 © d 2 

What is shown in [5 is that given a function con- 
structed in this way, one has: 



(out|C/|in> = ^Jl. 



(1) 



Here #o, #i 7 denote the number of solutions to the equa- 
tion f(x) = 0, f{x) = 1 respectively, and h denotes the 
number of Hadamard gates in the circuit. Note that 
#0 + #1 = 2 V where v is the number of variables in 
the function / once the input and output qubit values 
have been fixed. If there are q qubits in the circuit then 
v = h — q. 

There are several other points to note in terms of the 
construction of /. Firstly it will be convenient to assume 
that every variable goes through at most one Toffoli gate 
- this can be arranged by inserting double Hadamard (i.e. 
identity) gates where necessary. This should also be done 
at the final outputs to the quantum circuit. Doing so 
ensures that the function / has the following properties 
(i) it is (monotone) cubic (ii) every variable appears in 
at most one cubic clause and two quadratic clauses. 
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sum of the weighted cycle covers of the associated graph. 



Now counting solutions to a general GF(2) polynomial 
is a #P-complete problem [6 j. That is, it has the same 
complexity as computing the permanent of a matrix - 
the prototypical #P problem - as was famously proven 
by Valiant in 1979 [7]. So we know that in principle we 
can map between these problems and find some matrices 
M , Mi such that per(M ) = # , per (Mi) = #i and 
then 



(out | U | in) = 



per(Mp) -per(Mi) 
V2 h 



However; the actual mapping between these problems 
is not particularly simple or economical. In addition 
Valiant's construction of the matrix to count solutions 
of a satisfiability problem is also not particularly eco- 
nomical. 

The purpose of this article is to present a very sim- 
ple, direct and economical construction relating quan- 
tum computing to evaluating a matrix permanent, which 
is also considerably more efficient than following the pre- 
ceding route. Moreover; instead of expressing the solu- 
tion to the problem as the difference in two matrix per- 
manent s, we will construct a single matrix/graph G such 
that 



(outlC/lin) = E^. 



(2) 



The route to finding G uses some of the same tricks as 
in Valiant's proof. As this paper is intended to also be ac- 
cessible for physicists possibly unfamiliar with Valiant's 
result, we will try and make the presentation as self- 
contained as possible. 

Any n x n matrix can be considered the weighted ad- 
jacency matrix for a weighted graph on n vertices, where 
the weight on the edge between vertices i and j is simply 
the (i,j)'th element of the matrix. The permanent of a 
matrix, formally defined by 



per(M)= U M ^ 



with S n the symmetric group on n symbols, is then 
graphically equivalent to the sum total of the weighted 



cycle covers of the graph: A cycle in a graph is a closed 
path, a cycle cover is a set of cycles for which each ver- 
tex belongs to one and only one cycle. The weight of a 
cycle cover is the product of the weights on the edges in- 
volved in that particular cycle cover - so the permanent 
is the sum of all such weights. An example is provided 
in Fig. [2] A brief summary of how the permanent arises 
in some physical considerations can be found in [8]. 

Let us first give an overall view of the construction. We 
will be constructing a graph in such a way that the pres- 
ence/absence of one particular cycle in any given cycle 
cover corresponds to whether a particular boolean vari- 
able Xi associated with this cycle is or 1. We will use 
the convention that if the particular cycle is present in 
the cycle cover then this matches the variable assignment 
Xi = 0, if it is not then xi = 1. Not all cycles within the 
graph will correspond to variable assignments - the ones 
which do we term external cycles. In figures the "exter- 
nal edges" which can make up such cycles will be colored 
in blue (to aid the eye only - there is no mathematical 
difference between these edges and other edges in the 
graph). The overall graph will consist of some "graph 
gadgets" (small subgraphs) connected by external edges. 
An example is given in Fig. [3] Each of the gadgets cor- 
responds to a clause - in the figure we show only the 
vertices of the gadget which connect to external edges. 
The blue external edges form loops around two or three 
of the graph gadgets according to whether the variable 
appears in two or three clauses, and obviously they loop 
through a clause gadget with their corresponding part- 
ners of that clause looping through the other vertices of 
the gadget. 

Now as we compute the sum of the weighted cycle cov- 
ers of the graph (ie the permanent of the associated ma- 
trix) each cycle cover in the sum corresponds to a partic- 
ular assignment of values to the Boolean variables - i.e. 
it will have a particular set of external cycles traversed, 
setting those variables to a value 0. The graph gadgets 
will be designed so that if none of the external edges con- 
nected to that gadget are traversed - corresponding to all 
of the variables in that clause being equal to 1 - then the 
weight which that gadget contributes to the particular cy- 
cle cover is —1. In all other cases the weight contributed 
by that gadget will be +1. Recall that the weight of any 
given cycle cover is the product of the weights over all cy- 
cles in the cover. So for a fixed cycle cover (corresponding 
to a fixed assignment to the boolean variables) the total 
weight will be +1 or —1 according to whether an even or 
an odd number of clauses are satisfied by that particular 
assignment. Assuming w.l.o.g. an even number of clauses 
in total, this in turn means that the weight of the partic- 
ular cycle cover is +1 if f(x) = and — 1 if f(x) = 1. As 
we sum over all weighted cycle covers we automatically 
are calculating the difference in the number of solutions 
of f(x) = to f(x) = 1, which is precisely what we need 
by Eq. qamp. 

The inner working of graph gadgets which act in the 
desired manner are shown in Figs. ??. If in some cy- 
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FIG. 3: (Color online) A big-picture view of the construction. 
The external edges form loops through the graph gadgets, 
and each such cycle is associated with one particular Boolean 
variable x% . If the cycle is traversed in a particular cycle cover 
then that corresponds to setting that particular variable to 0. 
Conversely, if the particular cycle is not traversed then this 
corresponds to the associated variable having a value of 1. 
The graph gadgets have 2 or 3 vertices connecting to external 
edges according to whether they are gadgets for a quadratic or 
cubic clause. This graph would correspond to the polynomial 
X1X2X3 X1X5 X4X5X6 X5X6 X2X4 X6X7 X3X7 



cle cover no external edges are incident on the gadget 
then its contribution will be —1, which can be readily 
verified by computing the permanents of their adjacency 
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If one, two 



or three external edges are incident on the gadgets then 
the contribution to the cycle cover has weight +1, this is 
depicted in Fig. [4jc). 

There is one potential problem which has not been ad- 
dressed. What is to stop a particular cycle cover involv- 
ing only part of an external cycle corresponding to some 
given variable. Why, for example, do we not get screwed 
up by cycle covers which, say, enter at one vertex of the 
graph gadget but leave at a different one? A figurative 
picture of such an undesirable type of cycle is given in 
Fig.ga). 

The possibility of such problematic cycles is ruled out 
by the internal workings of the quadratic clause graph 
gadget. This is shown in Fig. [5jb) . Any cycle cover 
which enters the gadget along one external edge and tries 
to leave out via the external edge on the other side of 
the gadget has two possible paths for doing so. These 
paths pick up opposite signs, and so when summed over 
contribute to the total. The process is somewhat rem- 
iniscent of mach-zender interferometry! Note that we 
did not need to design the cubic graph gadget to have 
the same property. This is because in the formulation 
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FIG. 4: (a) The inner workings of the quadratic clause graph 
gadget. The weight on any edge is 1 unless otherwise indi- 
cated, (b) The inner workings of the cubic clause gadget. 
The weights (3i need only satisfy — —2, which can be 

achieved by setting /3i = = 03 = 1 if a graph with 

integer weights is desired, (c) The various ways in which the 
gadgets can be involved in a cycle cover with external edges. 
These are the cases when some of the associated boolean vari- 
ables in the clause are equal to and, as can be seen, these 
cases all contribute a weight of +1 to the cycle cover. When 
no external edges are incident on the gadget the weight it 
contributes must be -1, as discussed in the text. 




FIG. 5: (a) The sort of cycle covers we need to avoid cycles 
which only partially traverse an external cycle, as such setting 
the associated variable to in one clause and 1 in another, 
(b) The inner workings of the quadratic clause gadget which 
ensure that any external edge must exit by the same vertex 
it entered. The two depicted contributions to the cycle cover 
have opposite signs and cause the necessary cancelation. 
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FIG. 6: Putting everything together - how to draw the final 
graph G over the top of the associated circuit. Note that it is 
the variable on the target line of a Toffoli gate which is created 
by the Hadamard that acts after the Toffoli which actually is 
involved in the associated cubic clause. This graph would be 
computing the transition amplitude with |in) = 1 1 1 1 1) and 
|out) = 1 1100). 

we have chosen any variable appears in only one cubic 
clause, and it must then also appear in two quadratic 
clauses. The quadratic clause gadgets suffice to "force" 
an external edge which is incident into the cubic clause 
gadget to leave via the same vertex it entered. 

In terms of the basic construction the final thing to 
mention is that it is simple to force the values at the 



boundaries (the input/output to the circuit) to be 1 or 0. 
This is done either by simply not connecting any external 
edges into the associated gadget (setting the variable to 
1), or by forcing an external edge through the gadget by 
having that edge also loop through a vertex which has 
no "self loop" (setting the variable to 0). An example of 
this can be seen in Fig. [6] where the input qubits are all 
fixed to have value 1, and the top two qubits have value 
1 at the output while the bottom two qubits are set to 
the value of at the output [10 . 

The overall construction can be naturally laid out by 
drawing the graph directly on top of the circuit diagram. 
This is illustrated in Fig. [6] for the same circuit of Fig.[l] 

Note that the number of vertices in the graph G we 
associate to an given circuit is basically 3 times the num- 
ber of gates in the circuit. Let us denote this number of 
vertices as m. We have that 

If it were the case that || Jh/m II < 1 then the results of 

[9] imply there would exist an efficient classical algorithm 
to simulate this quantum circuit. 
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